Rotational Isotropy Breaking as Proof for Spin-polarized Cooper Pairs in the 

Topological Superconductor Cu x .Bi 2 Se 3 
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In a promising candidate of topological superconductors, Cu^I^Sea, we propose a way to ex- 
clusively determine the pairing symmetry. The proposal suggests that the angle dependence of the 
thermal conductivity in the basal afe-plane shows a distinct strong anisotropy only when the pairing 
symmetry is an odd-parity spin-polarized triplet below the superconducting transition temperature 
(2c)- Such striking isotropy breaking below T c is explicitly involved in Dirac formalism for super- 
conductors, in which the spin-orbit coupling is essential. We classify possible gap functions based 
on the Dirac formalism and clarify an origin of the isotropy breaking. 

PACS numbers: 74.20.Rp, 74.25.Op, 74.81.-g 



I. INTRODUCTION 

The discovery of topological insulator has received con- 
siderable attention because of its topologically protected 
nature of the surface gapless state. The non-trivial topol- 
ogy of the bulk insulator wavefunction shows up in edge 
boundaries together with closing the insulating gap 1 - - — . 
A typical topological insulator, Bi2Se3 is characterized by 
a non-trivial Zi topology on its bulk valence band, which 
brings about a helical edge state leading to spontaneous 
spin current. Theoretically, such a non-trivial feature can 
be also maintained by opening of superconducting gap. 
Presently, a quest for the so-called topological supercon- 
ductor is one of the most exciting issues in condensed 
matter physics. 

Very recently, owing to its close vicinity to Z% topologi- 
cal insulator Bi 2 Se3, Cu intercalated material, Cu K Bi 2 Se3 
has been regarded as a key compound to investigate the 
non-trivial topological superconductivity^ - — . Indeed, 
Cua;Bi2Se3 is a carrier doped compound, whose super- 
conducting transition occurs around 0.3K. Soon after the 
discovery, several groups observed zero-bias conductance 
peaks (ZBCP's) by using the point contact spectroscopy. 
Generally, ZBCP has been well-known as a signature of 
unconventional superconductivity such as d-wave or sign- 
reversing s-wave superconductivity^^. On the other 
hand, the origin of ZBCP observed in Cua;Bi2Se3 may 
be ascribed to topologically-protected gapless Majorana 
fermion at edges 1 ^— as a signature of topological su- 
perconductor. Then, clear evidence for such a protected 
edge state is now in great demand. 

So far, a tremendous number of studies have supported 
that the ZBCP is originated from Andreev bound-states 
formed at the edge boundary in unconventional super- 
conductors. The emergence of the bound states is deeply 
associated with the internal sign change in the unconven- 
tional Cooper pair. In the superconductor, Cua;Bi2Se3, 
Sasaki et al. theoretically demonstrated how the Ma- 
jorana fermion brings about ZBCP's- 2 . They examined 
four types of superconducting gap functions selected by 
the point-symmetry analysis and drew a theoretical re- 



mark that an odd-parity spin-triplet is the most-likely 
pairing symmetry to explain the observed ZBCP's. How- 
ever, they could not exclude the other paring symmetries, 
since any other three possible odd-parity pairing, one of 
whose gap is full and two of whose gaps have point- nodes, 
can induce ZBCP's with tiny variation of parameters in 
the original tight-binding model. Such confusing ambigu- 
ities reflect that it is very hard to identify the paring sym- 
metry only through the energy dependence of the density 
of states. In this paper, we therefore propose that angle 
dependence of thermal conductivity is a crucial probe to 
break the above controversy. 

In the history of the quest for superconducting pair- 
ing symmetry, the thermal conductivity has been fre- 
quently employed as a tool to identify the superconduct- 
ing gap structure. Its angle- dependence obtained by ro- 
tating the applied magnetic field allows to explore the 
variation in the gap amplitude implying the existence of 
gap node o 23 i 24 . On the other hand, in topological super- 
conductors, we do not need the application of the mag- 
netic field, since the expected Majorana bound-state at 
edges is protected by time-reversal symmetry. In this pa- 
per, we present that anisotropy of the thermal conduc- 
tivity is an exclusive proof to identify the paring state 
in the superconductor Cua;Bi2Se3. The odd-parity spin- 
polarized Cooper-pair breaks the horizontal-angle invari- 
ance of the quasi-particle eigen-states around T-point, 
while non-polarized ones preserve the rotational invari- 
ance. The consequence is demonstrated by numerical cal- 
culations on the quasi-particle thermal conductivity and 
supported by theoretical analysis on a low-energy Dirac 
formalism for the superconductivity. Our proposal does 
not at all depend on the material parameters of the tight- 
binding model for Cu x Bi2Se3, since the consequence is 
mathematically inherent in the formalism as long as its 
effective low-energy model is given by the Dirac-type one. 



II. MODEL AND METHOD 

We start with a model Hamiltonian on the topological 
insulator Bi2Se3 proposed by several groups, which in- 
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eludes two-orbital spin-orbit coupling in 4 x 4 matri x 12 ' 26 . 
In order to examine the edge state in the model with uni- 
form superconducting gap, we have the mean-field Hamil- 
tonian on discrete N z planes stacked along ^-direction (c- 
axis) based on the Bogoliubov-de Gennes formalism (see 
Fig. ED), 
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where e, is the 8-component annihilation operator at the 
i-th plane, and k x and k y denote the in-plane momentum. 
Then, we have 8x8 matrix Hamiltonian expressed as 
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( ^ij [kx-, ky) 



(2) 



where, A is 4 x 4 matrix whose elements are given as 
A^:, using the orbital l{m) and spin u(u') indices. The 
normal-state Hamiltonian ^{k x ,k y ) is written as 



{kx i ky) 
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Aij(kx,k y ) e i Ak x ,k y )so 
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where, so is the unit matrix in spin space, and ef^(k x , k y ) 
is expressed as ef-{k x ,k y ) = E±(k x ,k y )S i:j + (~Di + 
Bi)(S ij+ i + <%-i)/2. E± ee e(k x ,k y ) ± M(k x ,ky), 
where e(k x ,k y ) = 2D\ + D2rj{k Xl k y ) — /i, M(k x ,k y ) = 
M - 2B X - B 2V (k x ,ky) and f]{k x ,k y ) = (4/3) (3 - 
2 cos{k x ^f?>/2)) cos(fcj,/2 — cosk y ). The 2x2 matrix 
Aij(k Xl k y ) is given by 
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i?, whose elements 
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with the matrix 
pressed by [R a ] i:j = -iS ij+1 /2 + i<%_i/2 gen- 
erated by Fourier transformation of sink z and 
Af(k x ,k y ) ee (2/3)I 2 (v / 3sin(fc xV / 3/2)cos(fc y /2) + 
±i(cos(fc xV / 3/2) sin(^/2) + ain(k y ))). We set M = 0.28 
eV, ft = 0.5 eV, A x = 0.32 eV, A 2 = 4.1/a eV, 
B 1 = 0.216 eV, B 2 = 56.6/a 2 eV, D 1 = 0.024 eV, 
D 2 = 19.6/a 2 and a = 4.076 A as the material pa- 
rameters for Cu a ;Bi2Se3 12 . According to Refs^ i 12 ' 27 , 
we examine four different types of the superconducting 
gap function, Ai to A4, which cover possible all gap 
functions selected by the point-group symmetry analysis 
(see also Table [J). 

In order to obtain the angle dependence of the thermal- 
conductivity, we calculate the quasi-particle thermal con- 
ductivity tensor Kij(T) expressed a a 28 ' 29 



dE" 



(5) 



in which v™(k), E^, and t(E%), respectively, are quasi- 
particle velocity, energy, and relaxation time as a func- 
tion of k and the band index n, and / is Fermi distribu- 
tion function. For simplicity, we take r(Ek) = t. Then, 



we calculate the thermal conductivity tensor Eq. ([5]) by 
diagonalizing 8N Z x 8N Z matrix BdG Hamiltonian ©. 
It is noted that quasi-particles are localized around the 
top edge surface as shown in Fig. [T] We emphasize that 
the calculation result does not qualitatively change, 
even if one takes into account the other contributions to 
the thermal conductivity. This is because a key feature is 
whether the quasi-particle eigen-state, i.e., the electronic 
structure is horizontally angle dependent around T-point 
or not. 
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FIG. 1. (Color online) The schematic figure of the calculation 
target system, which is discretized along z-direction. Then, 
the quasi-particles are localized around the edge surface plane 
at z = 0. 



III. NUMERICAL RESULTS 

In order to examine whether the isotropy of K%j is 
preserved or not, we obtain the temperature depen- 
dence of k xx and K yy well below T c . Here, we drop 
the temperature-dependence of the pair amplitude as 
A(T) cx Ao = 0.05 eV for convenience of calculations. 
Then, the result is valid in the low temperature range as 
shown in Fig. [2j The thermal conductivity from the bulk 
body mainly arising from phonon is always isotropic free 
from the present argument. As shown in Fig. [51 the ther- 
mal conductivity is found to be isotropic and anisotropic 
reflecting the difference in the gap function. For the gap 



function A4 Q (Aj 2 



- A 12 - 



-A 21 - 



-A 21 . See, Ta- 



ble HI, one finds a clear anisotropy, while not for A 2 . 
This difference is explained by the horizontal-angle de- 
pendence of diagonalized (edge-state's) eigenvalue distri- 
bution as shown in Fig. [3] These results indicate that 
the gap function A^ a breaks the rotational isotropy in 
the original model without the gap function while other 
ones do not show any isotropy breaking. 



IV. ORIGIN OF THE ROTATIONAL ISOTROPY 
BREAKING: MASSIVE DIRAC EQUATIONS 

Now, let us theoretically pursue the origin of the above 
result. First, using a Nambu space not commonly used 
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FIG. 2. (Color online) The temperature dependence of the 
thermal conductivity carried by the edge bound states for the 
gap functions (a) A2 and (b) Ai a . N z = 64. 



in condensed matter theories, we introduce the massive 
Dirac Hamiltonian with superconducting pair function 
(see Appendix A.}^>2L22 



H = dr( $(r) Mr) 



h-(r) A~(r) 
A+(r) h+{r) 



with 



h± = M - t 1 ^ - j 2 id y - 7 3 z ± /i 7 °. 
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Here, A" = 7 °Ai 7 2 , A+ = 7 °(A-)V, and 7 l is_4 x 4 
Dirac gamma matrix in the Dirac representation. ip( x ) — 
V^(x) 7 °, V» c (a;) = ■0 < t 7 °, = C^ T , where C(= j 7 2 7 °) is a 
representative matrix of charge-conjugation 32 . The space 
is three-dimensional (3-D) x-y-z coordinate system. Uti- 
lizing 2x2 Pauli matrices dj in the orbital space and 
Sj in the spin space, gamma matrices are represented as 
7 ° = <j z ® 1, j l=1 ' 2 - 3 = ia y (g) ^ and 7 5 = a x ® 1, respec- 
tively, with the relation 7 5 = i 7 ° 7 1 7 2 7 3 . We note that 
parity of 7 ^° j s oc id since the transformation property 
is equivalent with that of vector. 

We classify the possible gap functions with gamma ma- 
trices in Table HI Reminding the symmetry of A(fe) = 



A(— fe), we have only six A matrices constrained by the 
fermion anti-symmetric relation — A = A T , 



AcxC,C 7 5 ,C 7 M 7 E 



(8) 



with = 0,1,2,3. According to Table [J (a list of A), 
the gap functions A~ introduced by Fu and Berg^21 are 
characterized by a scalar, a pseudo-scalar, and a four- 
vector: 



1 



(n = 2) 
(n = la) 
(else) 



(9) 



where, the Feynman slash fit is defined by X^7 Mq; a" 
and the gap function is characterized by the four-vector 
a. For example, the gap function A^ a is character- 
ized by — 5^2- The s-wave pairing gap function 
Ai Q in this representation has been well-studied in color 
superconductivity 33 . From the above compact expres- 
sions for the gap functions, one finds that the present 
Nambu representation using the charge-conjugation op- 
erator is quite convenient in the massive Dirac equations 
for the superconductivity. In terms of this representa- 
tion, we also show that gap functions represented by fic 
with 7 At ^°-components have point-nodes and know the 
positions of these point-nodes (see Appendix B). 

In the normal state, the Hamiltonian is rotationally in- 
variant in 2-D xy-space. On the other hand, in the super- 
conducting state, the gap function represented by fit with 
7 M ^°-components breaks the rotational symmetry, indi- 
cating non-zero spatial elements of four-vector a. Thus, 
the diagonalization together with A breaks the isotropy 
of quasi-particles. Of course, one easily finds, before the 
diagonalization step, that the gap function A2 maintains 
isotropy and full gap while A4 does not. Thus, as shown 
in Fig. O the distribution of eigenvalues of Eq. @ is 
horizontally angle-dependent only in the case of A4. 



V. DISCUSSION 

Finally, let us discuss the feasibility of the anisotropic 
thermal conductivity as a probe of pairing state. The 
first issue is quantitative predictability. In the present 
theoretical treatment, we neglect multi-orbital effects on 
quasi-particle scattering for simplicity and handle only 
orbital-diagonal one. Then, it means that quantitative 
predictability of the thermal conductivity is beyond the 
present scope. On the other hand, the anisotropy caused 
by the gap function A4 is still robust, since the quasi- 
velocities v x and v y become anisotropic below T c ow- 
ing to the coupling with A4. Clearly, such a contri- 
bution is a leading one in sufficiently low-temperature 
range T/Ao < 1 as shown in Fig. [2 where the bulk 
thermal conductivity irrelevant to the present mecha- 
nism is significantly reduced because of the supercon- 
ducting gap opening. The next is qualitative comparison 
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FIG. 3. (Left) The distributions of quasi-particle energies in in-plane momentum space with N z = 64 for the gap functions (a) 
A2 and (b) A4 a , respectively. (Right) The contour plots of the minimum eigenvalues regarded as the edge states for (a) A2 
and (b) A4 a , respectively. The contours are drawn at regular intervals (5 A — Ao/20). 
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TABLE I. The representation of the gap functions with gamma matrices. See text for A„. "P-Scalar" denotes a pseudo scalar 
whose parity is odd and "i-Polar" denotes a polar vector in i-direction in four-dimensional space. d v denotes a d-vector in the 
^-direction {[d v ]i = Ao<5j^). di and s% denote 2x2 Pauli matrices in the orbital space and the spin space, respectively. 



with other gap cases. In non-topological superconduc- 
tivity Ai a (A}| = -A|| = A^2 = -AJ2) the thermal 
conductivity is rather small and isotropic because of no 
edge bound state in contrast other gaps. For A 2 and 
A 3 (AjJ = Aj 2 . = -A^l = -A^), the gapless edge 
bound states might appear depending on the band pa- 
rameters, but the contribution to the thermal conductiv- 
ity is isotropic since these gap functions are symmetric 
around r-point in k x -k y plane. Thus, we point out that 
the isotropy breaking of the thermal conductivity is suffi- 
ciently an exclusive evidence of the spin-polarized Cooper 



pair (A 4 ). 

Moreover, our proposal does not depend on the ma- 
terial parameter values such as Mo, /1, A\, A2, B\, B2, 
Di, and D 2 - On the other hand, it should be noted 
that non-electronic disorder or other scattering contribu- 
tions might recover the original isotropy. However, the 
isotropy breaking mechanism is inherent in the topolog- 
ical superconductivity modeling. This idea is applicable 
for not only various topological superconductors but also 
superconductivity emerged in high-energy physics. 
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VI. CONCLUSION 



superconductors . 



In conclusion, we numerically calculated the thermal 
conductivity tensor k xx and K yy dominated by the edge 
bound states in the superconductor Cu^I^Ses. Con- 
sequently, we found that the rotational isotropy of the 
thermal conductivity is broken below T c for the spin- 
polarized gap A4. In order to explore the mechanism, 
we noticed that the Nambu representation expanded by 
charge-conjugation operator is convenient in handling 
massive Dirac equations with Cooper pair functions and 
newly classified the possible gap functions by using sim- 
ple mathematical tools as gamma matrices 7 and four- 
vector a. As a result of the classification, we found that 
the horizontal angle invariance of the quasi-particle eigen- 
values around T-point is broken for the spin-polarized 
pair A4, which is mathematically characterized by four- 
vector a lying on the basal xy-plane. We propose that 
the angle dependence of the thermal conductivity has an 
exclusive tool to identify the gap function in topological 
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Appendix A: Massive Dirac Hamiltonian 

We rewrite the BdG Hamiltonian to introduce the mas- 
sive Dirac Hamiltonian with superconducting pair func- 
tion ©. We start with the conventional BdG Hamilto- 
nian expressed as 



H = dr( V> + (r) VO) 3 



h(r)-fx 
At(r) 



A(r) 
-/i(r)*4 



(Al) 



r 



Here, the Hamiltonian in normal states is written as 

h(r) = Mq7° - l°l^d x - 7 °7 2 ^ - j°jHd z . (A2) 

We introduce the charge-conjugation operator, 

7P c {r)=Ci> T {r)=C 1 r{r), (A3) 

with the charge-conjugation matrix C = i7 2 7° and 
ijj(r) = i/j'(r)j°. The gamma matrices have the anti- 
commutation relation: 



yV + 7V = 2rf, 



(A4) 



with 77 = diag (1, — 1, — 1, — 1). We note that both the 
gamma matrices 7^ =1 > 2 ' 3 and C are the non-Hermitian 
matrices ( 7 M =1 . 2 .3t = _yi=i,2,3 and c t = _£). The ma- 
trix element ij} T {r)(— h*(r))ip*(r) is rewritten as 

i> T {r)(-h{ry)4,{ry = ^(r) 7 °tf (-fc(r)')C/^ c (r), 
= Mr)h(r)Mr)- (A5) 

Therefore, we obtain the massive Dirac Hamiltonian with 
the charge-conjugation operators ((5]). 



Appendix B: Positions of point- nodes 

In terms of the four vector a, one can easily find the 
positions of point-nodes. For example, the four vector 



for A^ a is characterized by a = (0, 0, 1, 0). In this case, 
the gap function has the rotational symmetry around y- 
axis. Therefore, the point-nodes can exist only on the k y - 
axis. Generalizing the above, we can show that gap func- 
tion represented by polar- vectors A ~ = Aoy =1,2,3 7 5 can 
have point-nodes. As we discussed in the above para- 
graph, the point-nodes can exist on the fc„-axis, since 
the gap function represented by a j/-polar vector has the 
rotational symmetry around j/-axis. Thus, we only con- 
sider the case of = 0. In this case, the normal- 
Hamiltonian in momentum space is written as 

h(k) = Mo7 + V7°Vy (Bl) 
The BdG equations are expressed as 

( 7 °Af + S^ Q k^ - fi)u + 7 V7 5 A « = Eu, (B2) 

7°7 I/ 7 5 A u + (7°M + 7%7 M + fj)v = Ev. (B3) 

If the point-nodes exist, the above equations have the 
zero-energy solutions (E = 0). By solving the BdG equa- 
tions, we obtain the relation, 



M 2 



(B4) 



On the other hand, the Fermi surface in normal states is 
determined by the following relation, 



K + K 



Mi 



(B5) 



Therefore, if the Ao is small, the point-nodes always exist 
near the normal-states Fermi surface. 
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